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We present a classification of energy flow variables for highly coUimated jets. Observables are 
constructed by taking moments of the energy flow and forming scalars of a suitable Lorentz subgroup. 
The jet shapes are naturally arranged in an expansion in both angular and energy resolution, allowing 
us to derive the natural observables for describing an A'^-particle jet. We classify the leading variables 
that characterize jets with up to 4 particles. We rediscover the familiar jet mass, angularities, and 
planar flow, which dominate the lowest order substructure variables. We also discover several new 
observables and we briefly discuss their physical interpretation. 



Analyzing the substructure of highly-boosted massive jets {pj > 500 GeV) has proved to be useful for 
distinguishing new physics signals from the QCD background. For instance, jets originating from boosted 
electrowcak gauge bosons [1 , tops [2 [3], Higgs [4|, and even new physics particles 016] are all interesting 
targets of searches conducted at the Tevatron and the LHC experiments, and it is therefore important to 
be able to distinguish them from QCD jets. For recent reviews on substructure techniques, experimental 
status and new physics searches see [7| and references therein. 

One way to characterize jet substructure is to consider observables which are functions of the energy 
flow within the jet, namely the energy distribution e{x) as measured by the detector, where x are the 
coordinates on the detector surface. Examples of such observables include the jet mass and the angularities 



where Ej is the jet energy, Ei is the energy of particle i, 6i is the angle relative to the jet axis, and a is a 
real parameter restricted to a < 2 to ensure IR safety. Another example is planar flow [9| (see also pi)]), 
which we discuss below in detail. 

While many observables have been defined, so far there has been no systematic way to classify them or 
to construct all observables at a given order in the detector resolution. The goal of the present work is to 
suggest such a classification for the case of narrow jets, as we now describe. 

An energy flow observable needs to meet several criteria. Ideally, it should be Lorentz invariant and 
IR-collinear (IRC) safe. In practice it is difficult to find observables which are Lorentz invariant (TU], the 
jet mass being an exception. 



I. INTRODUCTION 




(1) 



^ Here we use a different normalization than that of [8]. 
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We may however relax this requirement as fohows. For highly-boosted, narrow jets with fixed 4-momentum 
the cross section factorizes at leading order into the hard process and the jet function (in which we are 
interested). Since the jet momentum is an input parameter of the jet function, it is enough to consider 
observables which are invariant under the subgroup of Lorentz that doesn't change the jet momentum. 
For a massive jet, this little group is isomorphic to S0(3), as can be seen by boosting to the jet rest frame. 

Another property of a jet is that it is contained in a cone of radius R ^ 1 around its axis, and we should 
further restrict ourselves to Lorentz transformations that leave the cone invariant^. The subgroup that 
stabilizes both the jet momentum and the cone is SO (2), the rotations around the jet axis. 

In the present paper we focus on highly-collimated jets, working at leading order in the small cone size 
6i<R<^\. We will further assume that all the particles that make up the jet are massless. Under these 
assumptions, we propose a complete classification of energy flow observables which are IR safe and Lorentz 
invariant, up to Lorentz transformations that change the jet momentum and the cone. Furthermore the 
observables may be naturally arranged as an expansion in terms of the energy resolution and the cone 
size (or alternatively the energy and angular resolutions). 

In the narrow cone approximation the experimental calorimeter information of the jet is fully specified by 
the energy distribution on a surface perpendicular to the jet axis, e{x), with x = (xi, X2) corresponding to 
a set of coordinates on the two dimensional jet surface. To obtain the classification, we begin by describing 
a given energy distribution e{x) in terms of its moments. 



— j ^ X£"(x) Xi-^ ' ' ' Xi^ . 

Observables are then constructed by taking products of these moments, and forming S0(2) scalars by 
contracting their indices. Invariance under boosts along the jet axis can be achieved by adding a factor of 
Ej raised to some power. Since our moments form a straightforward generalization of angularities with 
a < 2, they are manifestly IR safe. 

Moments of higher rank are of higher order in the cone size i?, since x (x 9, the opening angle. Moreover, 
since / ~ e(a;), observables that are products of a larger number of moments are more sensitive to errors 
in energy measurement (we will be more precise about this point below). This allows us to arrange the 
observables in an expansion by energy resolution and cone size as will be shown in section |III[ At the 
lowest orders in this expansion we find the jet mass and the angularity with a = —2. At the next order 
we find the angularity with a = —4 and the planar fiow. Yet higher orders produce the whole tower of 
angularities with a = —2k, as well as many new observables. 

The paper is organized as follows. In the next section we review the definition of planar flow, showing 
how it naturally leads to the expansion of the energy distribution in terms of moments. In section III we 
define the moment expansion, find the first few observables, and precisely identify the small parameters in 
which we expand. In section [TV] we turn to a classification of the leading observables that characterize jets 
with up to 4 particles. In section [V] we present an alternative approach to the study of jet substructure, 
based on an expansion in terms of orthogonal functions. Finally, in section [Vl] we analyze one of the new 
observables discovered in the classification, developing techniques that may be applied to other observables 
as well. Section |VII| contains a few comments on the application of this formalism to the case of hadron 
colliders and our conclusions. 



^ Alternatively, we may describe the jet as its 4-momentum plus a detector surface of small area, perpendicular to the jet 
axis. The little group should leave the momentum and surface invariant. Moreover, even if we will refer to the jet cone 
throughout the paper, our conclusions are not restricted to jet cone algorithms. For us it is sufficient that there exist a 
number R < 1 such that each jet is contained in cones of radius R, which is also generically the case for sufficiently hard 
jets found with sequential recombination algorithms. Indeed, for narrow massive jets, it was shown that the angularity 
distributions are similar for a wide range of jet algorithms 
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II. MASS, PLANAR FLOW, AND THE SECOND MOMENT 

Planar flow is defined in terms of the matrix /.„ , with components'^ 

itEparticlcs it^particlcs 

where p± is the particle momentum in the detector plane, is the azimuthal angle, and fi{<j)) = 
cos(0) , /2((/') — sin(0) . The approximation is made under the narrow cone assumption. Planar flow 
pT| is then defined as 

4det/,„ 



(Tr/„)2 

One can easily verify that < Pf < 1, that it vanishes when the energy distribution lies on a line, and 
that it is maximal for isotropic distribution. 

Iiu also enters in the definition of the jet mass. From TOj — pj, pj — X^i-Pi it is easy to see that, in the 
narrow cone approximation, 

g « E -Tr/.. (2) 

i(i particles 

The narrow cone approximation ensures also that mj <^ Ej. 

Now, let Xk{9, (p) ~ OJk{(j)) denote a coordinate system on the detector plane. Then we may write I^j as 
a sum over detector cells, where we weigh each cell by its total collected energy E{x), namely 

cells 

itEparticlcs 

One can define a normalized version of this tensor, Iki — I*^ /Ej , where Ej is the jet energy. For the 
sake of clarity let us rewrite the above expression in integral form, 

Iki= J (fxe{x)xkXi , (3) 

where e{x) is the continuous, normalized energy distribution, given by 

e(a;) = — J2 E (^x^''^^ 5{x - x'^''^) . 

nGcells 

Note that 

d^xe{x) = l. (4) 



^ We use a different normalization for than that of 1111, without affecting the planar flow definition. 
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Equations ([3]) and Q are quite suggestive. The function s{x) encodes all the information about the jet 
structure and it is the object we want to characterize. Now, we have just seen that Iki is its second 
moment and it gives rise to two interesting physical observables — the jet mass and the planar flow. It 
is well known that all the information of e{x) is encoded in its moments. It is therefore plausible that 
expressing the energy distribution e{x) in terms of its moments would provide a natural way to derive 
observables with the desired properties. One can also expand the function e{x) in a set of orthogonal 
functions and characterize it in terms of its expansion coefficients. S0(2) invariance will be ensured 
by taking suitable coefRcient combinations that are singlets under rotation. The two approaches are 
complementary and can provide different insight on the properties of e{x). These are the ideas that we 
will explore in the remainder of this work. 



III. EXPANSION IN MOMENTS 



The n-th moment /„ of the energy distribution e{x) is defined by 

Iki,...,k„ = J d^xe{x)xki ■ ■ -Xk^ = ^ X! ^i^k]'"^kl- 

itzparticlcs 

The zeroth moment is Q. The first moment is the expectation value, or dipole. It is set to zero by the 
requirement that the total transverse momentum of the jet vanishes, a state which can be reached by 
rotating the jet. This fixes the origin of the detection plane, which in turn determines the jet axis. We 
have then 

/o = 1 , /i = . 

The first non-trivial moment is therefore I2, as expected. We are looking to define observables that are 
invariant under the little group SO (2), the Lorentz subgroup that doesn't change the jet momentum or 
the cone. S0(2) has two independent invariant tensors, 5ij and e^j. The only S0(2) scalar that is linear 
in I2 is the normalized version of ([2]), namely 

Next, consider a tensor product I2 <E) l2- There are three nontrivial scalars one may construct, 

^ii^jj : ^ij-^ij ' ^ij^kl-^ik^jl • 

Of these, only two are independent, since 

£i]<^kilikl]i = ^{liiljj - ifj) = 2det/. 

Also, the first scalar, luljj, factorizes in lower-rank scalars. We therefore find only one new observable, 
det(/), which is an un-normalized version of planar flow. 

Before proceeding with the expansion, let us clarify in what sense the moment expansion is an expansion in 
small parameters. Planar flow, composed of I2®l2i is apparently of a higher order than the mass squared, 
composed of a single l2- To quantify this statement, first note that energy distribution is constrained 
to lie within a small cone of radius i? ^ 1. Since /„ ~ ~ 0"', we have that /„ scales as i?". As an 
expansion in i?, planar flow is of order 4 while the jet mass squared is of order 2. 

Another small parameter we may consider is the angular resolution A0, which is ultimately limited by 
the calorimeter cell size. Including its effect gives /„ ^ (0 ± A6')" = 6*" ± nA6'6'"^^ + • • • . The error on 
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TABLE I: Moment products that correspond to given orders in energy resolution Ae and cone size 7?. Observables 
are constructed by contracting these products in various ways. 



the value of an observable due to the finite angular resolution is therefore proportional to the total rank 
of moments composing this observable. Planar flow has error ~ while the mass squared has error 

^ 2/S.d. Regardless of which small angular parameter we choose, measuring a moment of higher rank 
requires a more accurate detector. 

As for the energy, the only small parameter is the energy resolution Ae. A single moment is proportional 
to e(x), so for a given observable the error due to energy resolution will increase with the number of 
moments that make up the observable. Planar flow will tend to have a larger error than the mass. 



IV. THE FIRST FEW ORDERS 

We are now in a position to classify all the jet energy flow observables. 

Table |l] lists the outer products that appear in the first few orders of the expansion. Note that there are no 
observables with an odd power of R since we cannot fully contract an odd number of S0(2) indices. Each 
outer product may be contracted in several different ways, giving rise to different observables. However, 
algebraic relations will reduce the total number of independent contractions. 

On top of this expansion it is useful to consider jets that consist of a given number of particles (or detector 
cell towers). A jet of N particles is characterized by iN — 4 variables, corresponding to the number of 
particle momentum components, minus the jet axis, the jet energy, and the SO (2) angle associated to the 
overall rotation around the jet axis.* For given N we can then identify the lowest order, Lorentz-invariant 
observables that characterize such jets. Let us list the observables for jets with up to 4 particles. 

The two lowest order observables are 

lu^- y 

Ej ^ ' E^, ' 

i^particlcs 
iGparticlcs 

Therefore, the natural observables for describing a two-particle jet are the mass and the angularity t_2- 

To describe 3-particle jets we need 3 more observables in addition to the mass and r_2- As shown above, 
the outer product one independent contraction, eik^jilijlu ^ det/, the planar flow. Next we 



* This is at leading order. Soft phenomena, such as color connection with other jets in the event, will generally break this 
symmetry. 
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have /g which has one contraction, 

liijjkk — ^ ^ EiO^ = 32t_4, 

i^iparticlcs 

corresponding to another member of the angularity family with higher a. Generally, /2n with n > 1 has 
only one contraction, ^ J2i E^iO'^"', which corresponds to an angularity with a = 2(1 — n). 

At the next order we find /2^4, (-^3)^1 and Jg, with the following independent contractions: 

'■ ^ij^kllikljlmm 7 ^ij likljkll 
(-^3) ■ ^ij^kl-^ikm-Ijlni 7 ^ijk^ijk 
Is • hijjkkll 

The additional contractions lijiijkk, hjjlikk can be shown to be linearly dependent on these. Any 
combination of them can be chosen as the remaining observable for 3-particle jets. In lumping /§ with 
the rest we assumed for simplicity that energy and angular measurements have comparable weight in the 
small parameter expansion. Of course, if this assumption is not true, /§ may be preferred or disfavored 
in comparison with the other contractions. 

Note that an observable — such as planar flow — which includes one or more e symbols vanishes when 
all particles are on a line.^ Since two particles always lie on a line, these observables contribute at 
leading order to 3-particle jets. They can therefore be used to distinguish QCD jets, which favor 2-parton 
configurations, from e.g. top jets that favor 3-body decays. 

Finally, 4-particle jets are described by eight variables. So far we found nine leading order observables, 

hi 7 hijj 7 ^ij^kllikl jl 7 I-iijjkk 7 ^ij^kllikl jhnm 7 
^ij hk^jkll 7 ^ij^klhkml jlm 7 hjkhjk 7 hijjkkll • 

Of these, any leading eight can be chosen to describe the jet. 



V. EXPANSION IN ZERNIKE POLYNOMIALS 

One can expand e{x) in a series of orthogonal functions. Since e(r, ip) is defined on a disc of radius 
R, perhaps the most convenient expansion is in terms of the Zernike polynomials |12j . which form an 
orthogonal basis on the unit disc. They are defined by 

Z™(p,0) = i?™(p)cos(m</.), Z-"(p,0)-i?™(p)sin(m0), 0<p<l, 

where < m < n, n — m even, and -R™(p) are a set of polynomials of degree n respecting the orthogonality 
condition 

Jo 2n + 2 

The orthogonality among different m's follows trivially from the orthogonality of the Fourier modes. 
This set of functions is widely used in optics, in particular in the study of optical aberrations where the 
expansion coefficients have simple geometrical meaning. 



To see this, first rotate the fine configuration to he on the xi axis, and then note that the e tensor forces an X2 factor to 
appear. This factor vanishes wherever e{x) 7^ 0. 
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The expansion of the energy distribution is 

oo 



= ^ + XI [a«,mi?™(-^) cos(m(^) + a„ -mi?;^'^-^) sin(m(?!)) , (5) 

n=l 0<m<n, 
n—m even 

The conditions Iq ~ I and /i = fix ao,o = l/^i and ai^±i to vanish. One can further expand the 



moments defined in section fll in terms of the dn.m, finding that a moment of order r will be expressed 
as a linear combination of a„_rn's with n < r. Moreover, upon tracing over k of the r indices, m will be 
constrained to be < r — fc. An even/odd number of indices corresponds to m (and thus n) being even/odd. 

The lowest order invariants like the mass, the angularities with a = 2, 4 and planar flow have the following 
expression in terms of the Zernike coefficients an.m'- 

T-4 = l^-R ("6.0 + 7a4,o + 21a2,o + 35ao,o) , 

(i-pf)|f = ^i?M42 + 4-2)- 

In optics, the lowest order coefficients have been given names. In particular ai^±i is the tilt, a2_o is the 
defocus, ±2 are the 0° and 45°-astigmatism (respectively), 03,^1 is the coma, and 04,0 is the spherical 
aberration. The optical analogy may provide us with additional geometrical understanding of what is 
being probed by the jet shapes defined in the previous sections. 



VI. ANALYSIS OF ehh 



We now try to gain some intuition regarding one of the new observables found in section IV 
O = 2eijlikljkmm- Note that any observable which includes an odd number of e^- symbols is a pseudo- 
scalar with respect to parity on the detector plane. It therefore vanishes for any energy distribution that 
is symmetric under reflection through an arbitrary axis. 

Let I'j^ be the I4 moment traced once, namely Ikkij- Both I2 and I'^ are real, symmetric matrices, so we 
can write them in terms of the Pauli matrices, 

-'2 — ^2,0^ + ^2,1^ + -'2,2-^ , 
-'4 — ^4,0^ + ^4,1^ + -'4,2 Y ■ 

The observable can then be written as 

O = 2eijlikljkjnm, — 2Tr(/2e/4) = eijl2,il4:,j ■ 

Note that O doesn't depend on the cr*^ components (the trace). Treating /2,i,/4,j as d = 2 vectors, we see 
that 

= hxh, (6) 
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where the product is the cross product in d = 2 which produces a pseudo-scalar. 

To get further insight, let us compute the components of these vectors. We do this for a general moment 
^2k with all indices traced except two. It is expanded in Pauli matrices just like l2,I[ above, and we 
denote its corresponding d — 2 vector by T^f^. We write the result using polar coordinates, xi — rcos(0), 
X2 = r sin(^). 

Ik = Tr(/^,a3) = 7^,^11 - I^^,^ = ^ dr dc^e[r, </>) r^^+^ cos(2</>) , 

JQ JO 

Ilk = Tr(/^,ai) = 24,.i2 = dr dc^e[r, c^) r^^+^ sin(20) . 

Jo Jo 

Plug the Zernike expansion ^ of the energy distribution in l2k- The 4> integral picks out the m — ±2 
modes, and the radial integral picks up a linear combination of a2fc,±2, ■ • ■ , a2,±2, with + sign for the 
cosine and — sign for the sine. Specializing for the cases fc = 1, 2 we find: 

h^lR^i^'^A, = «4,2 + 5a2,2 y 

\a2~2 J 40 ^04, -2 + 5a2.-2 J 

The observable ^ can now be written as 

0=^R'det ""'A , (7) 

24U \a2,-2 04,-2 J 

where we have dropped the piece in I4 which is proportional to I2 because of the usual properties of the 
determinant. 

This shows that, in this expansion of £(2;), O depends only on the m — ±2 Fourier modes for the angular 
part and on two specific combinations of the even, n — 2,A Zernike modes. 

Finally, let us define the normalized observable 

On = 



I2 X h 
1/2 1 1/4 1 



which assumes values in the range [—1, 1]. This observable is maximal when I2 and I4 are orthogonal, and 
vanishes when I2, Ia are linearly dependent. This happens, for instance, when the energy distribution 
is invariant under reflection through some axis. Indeed, in that case we may rotate this axis to coincide 
with the Xi direction, following which a„__2 = since all the antisymmetric Fourier modes vanish. This 
agrees with the fact that this observable is a pseudo-scalar. 

Moreover, I2 and I'^ differ only by a different weighing of the inner/outer part of the jets (due to the 
additional in 14). A symmetric matrix defines a characteristic ellipsoid (in the case of I2 it is the 
ellipsoid of inertia in the detector plane) . A non-zero I2 x /4 determines by how much the principal axes 
of the ellipsoid rotate when we give extra weight to the outer portion of the jet. (In the language of optics 
it would correspond to measuring the relative orientation of the lowest and higher order astigmatism.) 
With this intuition it is also trivial to see that vanishes when all the particles lie at the same distance 
from the jet axis. 

We will now briefly study some features of the On distribution. While on an event-by-event basis On 
will take both positive and negative values, in most cases the distribution will be symmetric around 
zero, yielding (On) = due to its pseudoscalar nature. A necessary condition for (On) 7^ is that the 
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substructure of the jet is controlled by a parity violating interaction. However, the interesting study 
of whether and when (On) is non-zero is beyond the scope of this work; here we will merely present 
distributions for \On\ in a few characteristic kinematic scenarios. Since this observable vanishes for two- 
particle configurations, in Fig. [T] we present the distributions for three-particle final states. We include 
both the case of a pure one-to-three decay and the case of a three-body decay of particle A via two 
subsequent two-body decays, with an on-shell intermediate particle B (as would be in the case of top 
decays). In the latter case one can see the sharp edge, whose position is just determined by the ratio ^ 
mB/rriA- 




0.0 0.2 0.4 0.6 0.8 1.0 

|0„N|72X/4l/|72l|74l 



FIG. 1: The On normalized distribution for various values of iriB/rnA, the mother particle boost was fixed to 
EA/mA = 10. 



VII. DISCUSSION 



While the study presented in this paper is directly applicable to e'^e~ machines, let us briefly discuss 
its applicability to hadronic colliders, such as Tevatron and the LHC. Since the partonic center of mass 
system is unknown in a hadronic collision, one clearly has to have a definition for the moments which 
is invariant under boosts along the beam axis. At leading order in the narrow cone approximation this 
is trivially achieved by using A77 and A0 (differences in pseudorapidity and in azimuthal angle) for the 



We checked that it is insensitive to the values of R and of the jet energy Ej, if varied within reasonable ranges. 
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xi and X2 coordinates. In fact e(x) depends only on x and the ratios Ei/Ej. Since E = cosh 77 one 
finds that e{x) is also invariant at leading order in the narrow cone approximation. Corrections start at 
0(tanh77jA?7) and contaminate a moment of order r with moments of order k > r (due to the additional 
powers of A?;). However, being of higher order, they will be suppressed by powers of the cone size R, 
rendering the contamination small. 

Another difference with the e~^e~ case is in the structure of the intra-jet radiation. Indeed in terms of 
symmetries the presence of beam remnants which can be color-connected with the jet under study provide 
a specific direction, common to all events, of breaking of the SO (2) symmetry. While these effects have 
been widely studied in the literature in various contexts (see e.g O [131 [14] ) , it would still be interesting 
to see how they are represented in the language developed in this paper. While we will not embark here 
on a systematic study, we will sketch the analysis and outline a few results. 

One can relax the requirement of the S0(2) invariancc and still use the same formalism described in this 
paper to study these effects. In general there will be a set of directions corresponding to the beam axis, 
the directions of other jets in the event, etc. They identify two-dimensional unit vectors in the detector 
plane that can be used to contract the indices of the energy flow moment tensors. To lowest order and 
considering for simplicity the case of a single n one has 

Note also that in our counting of how many independent observables are present for an iV-particle system, 
S0(2) invariance eliminated one observable. Therefore, while the underlying event and soft radiation from 
color reconnection can contaminate all the jet shape observables presented in this paper, at leading order 
there is room for only one observable able to probe these effects via a systematic S0(2) non-invariance. 
This observable can be chosen to be any combination of the lowest order contractions such as the ones 
shown above. In particular it is worth noticing that lujU^ is closely related to t ■ n where t is the pull 
variable defined in |15| . the only differences being a different power of the cone size. 

To conclude, we have presented an order-by-order classification of jet shape variables for narrow jets with 
massless constituents. At the first few orders we encountered familiar variables — ^jet mass, angularity, 
and planar flow — and at higher orders we found new observables. Specifically, we proposed several new 
observables that may be used to characterize 3- and 4-particle jets. While the classification is complete, 
the formalism we introduced does not provide much insight into the geometric nature of the observables 
(beyond a few simple properties), and for that we need additional tools. Expanding the energy distribution 
in terms of Zernike polynomials seems especially suited for this purpose, at least for the set of observables 
analyzed here. It will be interesting to analyze other observables in this way, and in particular the 
techniques used for studying e/2/4 can be easily applied to other observables that involve the moments 
/2J.. A more thorough study of the properties of these new observables, and the identification of processes 
for which they provide discriminating power, is an important topic left to future study. 

The fact that the Zernike polynomials, which are commonly used in the field of optics, arise naturally in 
our jet shape analysis is obviously related to the fact that in both cases one is describing some distribution 
on a disc. One may hope however that the analogy between optics and jet substructure is deeper. If that 
is the case, further insight into the description of jet energy flow may be gained from the well-studied 
theory of aberrations. 
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